Mass ordering of differential elliptic flow and its violation for mesons 
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We simulate the dynamics of Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) 
with a hybrid model that treats the dense early quark-gluon plasma (QGP) stage macroscopically 
as an ideal fluid, but models the dilute late hadron resonance gas (HG) microscopically using a 
hadronic cascade. By comparing with a pure hydrodynamic approach we identify effects of hadronic 
viscosity on the transverse momentum spectra and differential elliptic flow «2(pt)- We investigate 
the dynamical origins of the observed mass-ordering of V2{pt) for identified hadrons, focusing on 
dissipative effects during the late hadronic stage. We find that, at RHIC energies, much of the 
finally observed mass-splitting is generated during the hadronic stage, due to build-up of additional 
radial fiow. The meson, having a small interaction cross section, does not fully participate in this 
additional flow. As a result, it violates the mass-ordering pattern for V2{pt) that is observed for 
other hadron species. We also show that the early decoupling of the meson from the hadronic 
rescattering dynamics leads to interesting and unambiguous features in the pT-dependence of the 
nuclear suppression factor Raa and of the (p/p ratio. 
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I. INTRODUCTION 

A presently hotly debated question is whether the 
quark-gluon plasma (QGP) created in Au-|-Au collisions 
at the Relativistic Heavy Ion Collider (RHIC) [l| repre- 
sents a "perfect liquid" 0, S, i, @, i.e. a fluid whose 
shear viscosity to entropy ratio rj/s is at or close to the 
conjectured [y] lower bound 7 = 1^. A key observable 
in this context is the elliptic flow V2 of hadrons emit- 
ted anisotropically in non-central collisions 0. At the 
highest RHIC energy of = 200 A GeV, the observed 
V2 values near midrapidity (jTyj^l), for not too large 
impact parameters (6^7 fm) and transverse momenta 
(pt ^ l.SGeV/c), agree well with predictions from ideal 
fluid dynamics [2] (i.e. assuming zero viscosity), includ- 
ing d, Q the predicted dependence of V2 on the trans- 
verse momentum pT and hadron masses [l0| . Evidence 
for non-zero shear viscosity in the collision fireball is ob- 
tained from deviations from ideal fluid dynamical behav- 
ior. This is manifested in the experimental data via a 
gradual break-down of the ideal fluid description when 
collisions are studied at larger impact parameters and at 
lower energies [TTj or when measurements are made away 
from midrapidity [H, [H, [13. In previous work [15[ we 
have shown that a large (and possibly the dominant) frac- 
tion of these deviations from ideal hydrodynamics is due 
to "late viscosity", caused by dissipative effects in the 
dilute hadronic rescattering stage that stretches between 
hadronization of the QGP and final kinetic freeze-out. 
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The question whether there is also non-negligible "early 
viscosity" in the dense QGP phase remains open. An 
answer to this question requires a proper viscous hydro- 
dynamic treatment of the QGP fluid which is presently 
being pursued vigorously [J, H, . It also depends on 
still unknown details of the initial conditions in heavy- 
ion collisions, in particular the initial spatial eccentricity 
of the fireball 0, [13, [H, [ll ■ 

In this paper we report additional results from the hy- 
brid model study presented in Ref. [Tsj . focusing our at- 
tention on a detailed investigation of dissipative effects 
during the late hadronic rescattering stage. The early 
QGP stage, including its hadronization, is described by 
ideal fluid dynamics. Speciflcally, we address the ques- 
tions of (i) how radial and elliptic flow evolve during 
the hadronic stage when it is described by a realistic 
hadronic rescattering cascade, rather than by an ideal 
fluid; (ii) how these differences affect the shapes of the 
finally observed hadronic p^-spectra and their differen- 
tial elliptic flow V2{pt)] and (iii) whether the differences 
between ideal fluid and realistic kinetic behavior during 
the late hadronic stage are similar for all hadronic species, 
or whether different magnitudes of their scattering cross 
sections translate into measurably different characteris- 
tics of their observed spectra and elliptic flow. 

The paper is organized as follows: For completeness, 
we begin in Section |TT] with a short review of the hy- 
brid model jl^l employed in this study. Our results are 
presented in Section IIIIl in three subsections organized 
along the questions raised in the preceding paragraph. 
We close the paper by presenting our conclusions and 
some perspectives in Section ITVl 
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II. THE MODEL 

Our study is based on a hybrid model which combines 
an ideal fluid dynamical description of the QGP stage 
with a realistic kinetic simulation of the hadronic stage 
. Relativistic hydrodynamics is the most relevant for- 
malism to understand the bulk and transport properties 
of the QGP since it directly connects the collective flow 
pattern developed during the QGP stage with its equa- 
tion of state (EOS). It is based on the key assumption of 
local thermalization. Since this assumption breaks down 
during both the very anisotropic initial matter formation 
stage and the dilute late hadronic rescattering stage, the 
hydrodynamic model can be applied only during the in- 
termediate period, between initial thermalization after 
about Ifm/c 0] and the completion of the quark-hadron 
transition which, in central Au-|-Au collisions at RHIC 
energies, happens after about lOfm/c. 

In absence of a non-equilibrium dynamical model for 
the early pre-equilibrium stage of the collision, its output 
is replaced by a set of initial conditions for the hydrody- 
namic evolution which are tuned to experimental mea- 
surements of the final state in central (6 = 0) collisions 
To describe the breakdown of the hydrodynamical 
model during the late hadronic stage, due to expansion 
and dilution of the matter, one has two options: One can 
either impose a sudden transition from thermalized mat- 
ter to non-interacting, free-streaming hadrons through 
the Cooper- Frye prescription [23|, imposed at a suitable 
value of the decoupling temperature Tdec or decoupling 
energy density Cdec 01j make a transition from the 
macroscopic hydrodynamic description to a microscopic 
kinetic description at a suitable value for the switching 
temperature Tg-w where both descriptions are simultane- 
ously valid [il, [m, m, m, [111, letting the subsequent 
kinetic decoupling play itself out automatically by follow- 
ing the microscopic evolution until all interactions have 
ceased. We here use both approaches alternatively, in 
order to isolate effects that are specifically caused by dis- 
sipative effects in the hadron rescattering cascade. 

A. Ideal Hydrodynamics 

For the space-time evolution of the perfect QGP fluid 
we solve numerically the equations of motion of ideal fluid 
dynamics, for a given initial state, in three spatial dimen- 
sions and in time ((3+l)-d ideal hydrodynamics) [isj : 

a^T^"^ = 0, (1) 

Tf"" = {e + p)u''u'' - pg'"' . (2) 

Here e, p, and u'^ are energy density, pressure, and 
four-velocity of the fluid, respectively. Due to its small- 
ness at collider energies, we neglect the net baryon den- 
sity. As an algorithm to solve the above equations we 
choose the Piecewise Parabohc Method (PPM) It 
is known to be a very robust scheme for solving non- 



relativistic gas dynamics including shock wave forma- 
tion and has been employed in many fields. We first 
applied it in [26| to solve Eulerian hydrodynamics for 
relativistic heavy- ion collisions, Eqs. ([T]) and Use 
of this algorithm enables us to describe the space-time 
evolution of relativistic fluids accurately even if the mat- 
ter passes through a first-order phase transition. The 
PPM is a higher order extension of the piecewise linear 
method employed, for example, in the rHLLE algorithm 
[13. We solve Eq. ^ in (r, x, y, f]s) coordinates [Tsj 
where r = \/t^—z'^) and 77^ = ^ ln[(t+z)/(t— z)] are lon- 
gitudinal proper time and space-time rapidity, respec- 
tively, adequate for the description of collisions at ultra- 
relativistic energies. The grid sizes are Ar = 0.3fm/c, 
Ax = Ay = 0.3 fm, and Ar/s=0.3. We have checked the 
grid size dependence of our final results and observed suf- 
ficient convergence with the given choice of grid param- 
eters, as long as smooth initial conditions such as those 
discussed below are used. 



B. Equation of State 

For the high temperature [T > = 170 MeV) QGP 
phase we use the EOS of massless non-interacting parton 
gas (u, d, s quarks and gluons) with a bag pressure B: 

p=i(e-4S) (3) 

The bag constant is tuned to = 247. 19 MeV to ensure 
a first order phase transition to a hadron resonance gas 
at critical temperature Tc = 170MeV. The hadron reso- 
nance gas model at T < includes all hadrons up to 
the mass of the A(1232) resonance. Systematic studies 
with various models of the EOS including a more realistic 
cross-over one will be discussed elsewhere. 

For a meaningful discrimination between the ideal fluid 
and hadron cascade descriptions of the hadron phase, 
and a realistic direct comparison of hydrodynamic results 
with experimental data, our hadron resonance gas EOS 
implements chemical freeze-out at Tchom = Tc = 170 MeV, 
as observed in RHIC collisions [2^. This is achieved by 
introducing appropriate temperature-dependent chemi- 
cal potentials Hi(T) for all hadronic species i in such 
a way that their numbers Ni including all decay con- 
tributions from higher- lying resonances, Ni — Nj + 

So h R ^ixNg . are conserved during the evolution |14l . 
[El, [HI, m. [Here N, is the number of the i- 
th hadron, and bu—^ix is the effective branching ratio (a 
product of branching ratio and degeneracy) of a decay 
process R ^ i + X.] In this "PCE model" [l4| only 
strongly interacting resonances with large decay widths 
(whose decays do not alter Ni) remain chemically equili- 
brated below the chemical freeze-out temperature. 

The hadronic chemical composition described by hy- 
drodynamics using the PCE model EOS is roughly con- 
sistent with that of the hadronic cascade models, as long 
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as the latter are initialized at Tgw with thermal and chem- 
ical equilibrium distributions [3^ . This is crucial for a 
meaningful comparison between hydrodynamic and ki- 
netic descriptions of hadronic matter since the chemi- 
cal composition of the hadron resonance gas has a sig- 
nificant influence on the hydrodynamic evolution of the 
hadronic transverse momentum spectra While the 
non-equilibrium hadronic chemical potentials UiiT) do 
not affect the EOS p{e) of the hadronic phase and 
thus lead to almost identical evolution of radial flow and 
total momentum anisotropy as for a chemically equili- 
brated hadron gas, they significantly alter the relation- 
ship between energy density and temperature, leading to 
cooler temperatures and hence to steeper transverse mo- 
mentum sgectra at identical kinetic decoupling energy 
densities . This effect is seen most dramatically in the 
time-dependence of the mean transverse momentum for 
pions [3|: {pt)tt decreases with proper time after chem- 
ical freczc-out whereas with continued hadronic chemi- 
cal equilibrium it would increase with time. Clear con- 
clusions about hadronic dissipative effects on the shapes 
of the transverse momentum spectra can therefore only 
be drawn from a comparison with hydrodynamic mod- 
els that implement chemical and kinetic freeze-out sepa- 
rately. 



C. Initial Conditions 

Contrary to Ref. [lB| where we studied both Glauber 
model and Color Glass Condensate (CGC) type initial 
conditions, for the comparative study presented here we 
concentrate on the Glauber model, suitably generalized 
to account for the longitudinal structure of particle multi- 
plicity (isl . [36| . We assume an initial entropy distribution 
of massless partons according to 
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where x± = (x, y) is the position perpendicular to the 
beam axis, C is a normalization factor, the "soft fraction" 
a is explained below, the parameter Yb is the beam ra- 
pidity, and is a suitable parametrization of the shape 
of rapidity distribution in pp collisions: 



r{r,s) = exp 
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number of wounded nucleons in each of the two nuclei 
and the number of binary nucleon-nucleon collisions, re- 
spectively. These are calculated from the Glauber model 
nuclear thickness function Ta,b{x±) [38| : 
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Here (j'§j^ ~ 42 mb is the inelastic nucleon-nucleon cross 

' 2/^1 (where b is the im- 



section, and r± ~ [(^±^6 
pact parameter). 

The soft/hard fraction a = 0.85 was adjusted to re- 
produce the measured centrality dependence [soj of the 
charged hadron multiplicity at midrapidity. At ?7s = 0, 
Eq. Q reduces to ^^/J^^ cx 



(1— a)ncoii| where n; 



dN 



[40| ; this parameterization is 



equivalent to the one used in [4l|, oc ("-part+^-part) + 
a^ncoii] , with x = From Eq. we can compute the 
entropy density at the initial time tq =0.6fm/c of the 
hydrodynamic evolution, s{tq, a;_L, rjs) ^ 



dS 



which 



TQdr]sd'^x ± ' 

provides the initial energy density and pressure distribu- 
tions through the tabulated EOS described above. 

Glauber model initial conditions have a long tradition 
for hydrodynamic simulations of heavy-ion collisions. In 
our previous study [isj we showed that with such ini- 
tial conditions "late viscosity" effects during the dilute 
hadronic rescattering stage are sufficient to explain all 
observed deviations of elliptic flow measurements from 
ideal fluid dynamical predictions. No significant addi- 
tional viscous effects during the early QGP stage were 
necessary. We also noted, however, that this conclu- 
sion depends crucially on this particular choice of ini- 
tial conditions, specifically the initial source eccentric- 
ity predicted by the Glauber model. The good agree- 
ment between theory and experiment disappears when 
one instead calculates the initial conditions from the KLN 
model [13, m, [H, [U, m, m, which is based on CGC 
ideas and, for the same impact parameter, produces al- 
most 30% larger source eccentricities. If Nature gives 
preference to such more eccentric initial conditions, ad- 
ditional viscous effects and/or a softer EOS for the QGP 
stag e may be needed to reproduce the experimental data 
[Ud^l. Here, we will not pursue this line of thought 
any further, but focus on the case of Glauber model ini- 
tial conditions and the specific modifications of hadron 
spectra and flow caused by "late hadronic viscosity" . 



We study Au+ Au coUisions at y/s ~ 200 A GeV and use 
C = 24:, A?7 =1.3, and cr,, = 2.1, so chosen as to reproduce 
the charged hadron pseudorapidity distributions mea- 
sured in these collisions H^]. N^^^ and iVcoU are the 



D. Hadronic Cascade Model 

In our hybrid model simulations we switch from ideal 
hydrodynamics to a hadronic cascade model at the 
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switching temperature Tgw = 169 MeV. The subsequent 
hadronic rescattering cascade is modeled by JAM [4^ . 
initiahzed with hadrons distributed according to the hy- 
drodynamic model output, calculated with the Cooper- 
Frye formula [2^ along the Tgw = 169 MeV hypersurface 
rejecting inward-going particles. We have checked [lB| 
that switching from an ideal hydrodynamic to a hybrid 
model description does not entail a major readjustment 
of initial conditions: Keeping the same initial conditions 
and hard/soft fraction a as previously determined within 
a purely hydrodynamic approach (see [3, fill for a de- 
tailed discussion of that procedure) we find [la| that the 
ccntrality dependence of dNch/drj at midrapidity remains 
consistent with the experimental data even if we switch 
below Tsv, to the hadronic cascade. Effects on the hadron 
spectra and elliptic flow are significant, however, and will 
be discussed in the next section. 



As customary in hadronic cascade models [45|, |46|, l47[ , 
JAM implements experimental hadronic scattering cross 
section data where available and uses the additive quark 
model where data do not exist, assuming the following 
formula for the total cross section: 



Ctot 



3 3 



4°^^? f 1 - 0.4^ ) ( 1 - 0.4:^ ) . (9) 



n2 



Here is the total nuclcon-nuclcon cross section, nt is 
the number of constituent quarks in a hadron, and Usi is 
the number of strange quarks in a hadron. For hadrons 
composed entirely of strange quarks, such as = (ss) and 
n = {sss), the cross sections become very small, due to 
the suppression factors in brackets in Eq. Only when 
we calculate spectra for (j) mesons in Sec. lIII Cl the decay 
channels for (f> mesons are switched off in the hadronic 
cascade calculations. Since the life time of (j) mesons 
46 fm/c) is longer than the typical life time of the system 
10-20 fm/c), and the number of (j) mesons is small 
compared to pions, kaons, and nucleons, this prescription 
is not expected to affect the bulk space-time evolution 
during the hadronic stage. 



III. RESULTS 



In Ref. [1^ we investigated the effect of hadronic dis- 
sipation on elliptic flow and found that it significantly 
suppresses the p^-integrated V2 at forward and backward 
rapidity and in peripheral collisions. In the following we 
explore the origins of this finding in more detail, by in- 
vestigating hadronic dissipativc effects on hadron spectra 
and differential elliptic flow V2{pt)- We finally explore 
specifically the spectra and elliptic flow of mesons as 
an example of a hadron that is only weakly coupled to 
the rest of the expanding hadronic fireball. 



A. Hadronic dissipative effects on spectra and 
elliptic flow 

In this subsection, we compare results from the hy- 
brid model with the ones from ideal hydrodynamics. In 
ideal hydrodynamic calculations it is assumed that even 
the late hadron resonance gas phase is characterized by 
essentially vanishing mean free paths and thus behaves 
as a perfect fluid, all the way down to kinetic decou- 
pling of the hadron momenta at Tth^lOOMeV. (This 
value is obtained by a simultaneous flt of the pion and 
proton spectra in central collisions which allows to sepa- 
rate the effects of radial flow and thermal motion at ki- 
netic freeze-out 0.) As discussed, chemical freeze-out is 
implemented at Tchom = 170 MeV by using an EOS with 
non-equilibrium chemical potentials which hold the sta- 
ble particle yields constant (and close to the ones in the 
cascade model approach) during the hydrodynamic evo- 
lution of the hadronic phase. The key difference between 
the hydrodynamic and hybrid model approaches is, thus, 
the finite mean free path for momentum-changing colli- 
sions in the hadronic cascade. 

In Fig. [1] pT-spectra for protons and pions arc shown 
for both, the hybrid model and the ideal hydrodynamic 
approach. For comparison, we also plot the pT-spectra 
without hadronic rescattering, obtained by setting all 
cross sections to zero in the hadron cascade or by set- 
ting Tth ~ Tsw ~ 169 MeV in the hydrodynamic approach 
(both procedures give the same spectra, by construction) . 
Note here that we include contributions from all reso- 
nances (except for weak decays unless explicitly noted 
otherwise) in ideal hydrodynamic and hybrid-model re- 
sults. One sees that hadronic rescattering in the JAM 
cascade pushes the protons to higher pT in exactly the 
same way as the growing radial fiow does in the hydro- 
dynamic approach, if one chooses for the latter a kinetic 
decoupling temperature of Tth = 100 MeV. The reason- 
able fit of the measured proton p^-spectra 
Pt ^ l-5GcV /c by the hydrodynamic model 
thus persists in the hybrid model approach (see Fig. [3] in 
the following subsection). 

The lack of visible dissipative effects on the proton 
spectra is probably an artifact caused by a judicial choice 
of the kinetic freeze-out temperature Tth = 100 MeV in 
the hydrodynamic approach, which was driven by the 
wish to reproduce the measured proton spectra with this 
model. This accident does not repeat itself for the pi- 
ons, shown in Fig.[IJb). For pions, the px-spectrum be- 
comes slightly steeper when evolved hydrodynamically 
(the steepening effects due to cooling are not quite com- 
pensated by the increasing radial fiow) whereas it gets 
hardened by hadronic rescattering effects in the hybrid 
approach. 

This pattern is consistent with theoretical expecta- 
tions: In the ideal fluid approach, pdV work in the longi- 
tudinal direction reduces the transverse energy per unit 
rapidity [4^ . [sot . Since pions dominate the medium but 
their number is fixed after chemical freeze-out, this leads 
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FIG. 1: (Color online) pr spectra with and without hadronic 
rescattering for (a) protons and (b) pions at midrapidity for 
Au+Au collisions at b — 2 fm, compared with results from 
ideal hydrodynamics decoupling at Tth = 100 MeV. 



to a decrease of the average px per pion Q , explaining the 
steeper pion spectrum from ideal hydrodynamics. (This 
argument is not quantitative since it neglects the shifting 
balance of transverse energy carried by pions and heavier 
particles such as protons which are more strongly affected 
by the developing radial flow Q. Also note that it does 
not remain true if a chemical equilibrium EOS is used 
in the hadronic phase where the pion number decreases 
with temperature and the average transverse energy per 
pion thus increases Q.) In contrast to the ideal fluid, the 
hadron gas in the JAM cascade is highly viscous. Shear 
viscosity is known to reduce the longitudinal and increase 
the transverse pressure [l^, reducing the loss of trans- 
verse energy due to longitudinal pdV work and increasing 
the transverse flow due to larger transverse pressure gra- 
dients [l6j . In addition, there arc viscous corrections to 
the (flow-boosted) thermal equilibrium form of the dis- 
tribution function at kinetic frcczc-out which lead to an 
additional viscous distortion of the p^-spectrum which 
actually increases with jsij . For Bjorken expansion 
of a homogeneous cylinder this distortion can be written 



analytically as [5l| 
dN 



prdpq 



dNp 

ATfT'^^^ J prdpT 



(10) 



where is the spectrum calculated from a boosted 

thermal equilibrium distribution along the decoupling 
surface at freeze-out time Tf and temperature T, and the 
expression in brackets preceding it is the p^-dependent 
viscous correction, parametrized by the sound attenua- 
tion length = I ^ (where rj is the shear viscosity) . 

The viscous flattening of the pion spectrum relative to 
the pure hydrodynamic approach seen in Fig. [ijb) re- 
ceives contributions from both factors in Eq. pH)) : ^'^'^"^ 
is flattened by the larger transverse flow generated by 
the viscously increased transverse pressure, and addi- 
tional flattening comes from the factor in brackets, due 
to a non-zero value for F^ in a viscous fluid. Wc don't 
know which of the two effects is larger; we only note that 
the pion spectrum from the hybrid model can be fitted 
very well by simply multiplying the hydrodynamic model 
spectrum with the factor in brackets in Eq. ([TU|) . tak- 
ing r = Tth = 100 MeV and adjusting F,/t/=0.01. How 
meaningful such a fit is (given that the form (jlO|) makes 
unrealistic assumptions about the fireball expansion) re- 
mains to be seen when realistic viscous hydrodynamic 
studies become available. 
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FIG. 2: (Color online) V2{pt) for pions and protons in I77I < 1.3 
at b = 7.2 fm. Results for pions (solid) and protons (dotted) 
from ideal hydrodynamics with Tth = 100 MeV are compared 
with the ones for pions (dashed) and protons (dash-dotted) 
from the hybrid model. 

While these considerations provide a qualitative expla- 
nation for the harder pion p^-spectrum from the JAM 
cascade compared to ideal hydrodynamics, the same ar- 
guments should also hold for protons where no such ef- 
fects are seen in Fig. [T{a). As already stated, this is 
presumably a consequence of an accidental cancellation 
of delicate thermal and flow effects with viscous correc- 
tions for our specific choice of Tth in the hydrodynamic 
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FIG. 3: (Color online) Centrality dependences of the pT spectra for (a) pions, (b) kaons, and (c) protons obtained from our 
hydro+cascade hybrid model, compared with data from the PHENIX Collaboration [43| for 200 A GeV Au+Au collisions. 
Impact parameters are (from top to bottom) h = 2.0, 7.2, and 9.7 fm, corresponding to the 0-5%, 20-30%, and 30-40% centrality 
ranges, respectively. 



model. Again, a full understanding of these results may 
require comparison with a viscous hydrodynamic treat- 
ment 

Figure [5] shows the pt dependence of V2 for pions and 
protons in semi-central Au+Au collisions {h = 7.2 fm) at 
midrapidity (|77|<1.3), comparing results from the hy- 
brid model with ideal hydrodynamics. Whereas, after 
an initial quadratic rise which extends over a larger nr- 
range for the heavier protons than the lighter pions [lO[ , 
the differential elliptic flow V2{pt) from ideal hydrody- 
namics increases almost linearly with pT, this increase is 
tempered in the results from the hadronic cascade. The 
differences between the two models is seen to grow with 
increasing pT- Again, this is qualitatively just as ex- 
pected from shear viscous effects 0, i, [Hf. Obviously, 
the different transport properties of the hadronic matter 
in JAM and in hydrodynamics are seen more clearly in 
the differential elliptic flow V2{pt) than in the pt spectra. 

B. Spectra and elliptic flow for tt, K, and p 

In this subsection, wc compare our results from the hy- 
brid model with experimental data for identified hadrons. 
In Fig. [21 transverse momentum spectra for pions, kaons, 
and protons from the hybrid model are compared with 
PHENIX data [i^ , for three impact parameters (central- 
ity classes) as shown in the figure. (The impact param- 
eters are adjusted to give the correct average number of 
participants for each centrality class, as quoted in [48|.) 
In all cases, the experimental data are reasonably well 
reproduced by the hybrid model for low transverse mo- 
menta to Pt ^ 1.5-2.0 GeV/c. Additional components 
(such as thermal quark recombination and jet fragmen- 
tation, including energy loss of fast partons in the fireball 



medium) would be required to reproduce the data above 
Pt ^ 1.5 GeV/c. It should be emphasized that, unlike in 
the purely hydrodynamic approach where the pT slope is 
controlled by the choice of kinetic freeze-out temperature 
and the correct hadron yields are ensured by appropriate 
choice of non-equilibrium hadron chemical potentials (see 
Sec. Ill Bp . the hybrid model has no adjustable parame- 
ters to reproduce both slope and normalization of the 
transverse momentum spectra. Hadronic cascade pro- 
cesses automatically describe both chemical and kinetic 
freeze-out. 

In Figure |4l we compare the pt dependence of V2 for 
pions, kaons, and protons with the STAR data for f2{2} 
[54| . for four centrality classes. For the 0-5% centrality 
class we show only pions since the quality of the kaon 
and proton data at this centrality is insufficient for a 
meaningful comparison with theory. The hybrid model 
correctly describes the mass ordering of the differential el- 
liptic flow, V2{pt) > V2{pt) > V2{pt), as seen in the data 
within the low-px region covered by the figure. Quanti- 
tatively, it provides a reasonable description up to 50% 
centrality, except for the most central collisions: Our re- 
sult for pions at 6 = 2.0fm is significantly smaller than 
the data. This can be attributed to the absence of eccen- 
tricity fluctuations in our model calculations [l^, [11] . 

To better understand the origin of the mass ordering 
in V2{pt), we compare in Fig.[5l for a selected impact pa- 
rameter of 5 = 7.2 fm, the above hybrid model result with 
a calculation where all hadronic rescattering is turned 
off, allowing only for decay of the unstable hadron reso- 
nances. Whereas just after hadronization the differential 
elliptic flow V2{pt) for pions and protons looks very simi- 
lar, the mass splitting gets strongly enhanced by hadronic 
rescattering. The smallness of the pion-proton mass split- 
ting at Tgw is partially accidental, because the splitting 
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FIG. 4: (Color online) Transverse momentum dependence 
of the elliptic flow coeflicient 112 for pions (dotted blue), 
kaons (solid red), and protons (dashed green) from the hybrid 
model, compared with STAR data for V2{2} from 200 AGeV 
Au+Au collisions, in four centrality classes [53 |. 



caused by the radial flow already established during the 
hydrodynamic QGP phase is significantly decreased 
by the effect of resonance decays which reduces the pion 
elliptic flow 'uJ(pt) by about 15% [l^,!!!]. Hadronic evo- 
lution below Tsw steepens the slope of V2(j)t) for pions 
[3], due to the generation of additional (integrated) V2 
and the reduction of their mean transverse momentum 
{Pt}tt Hi- (Note that for pions the slope of V2{pt) can 
be simply approximated as dv2{pT) / dpx ~ V2/ (pt) 0-) 
For heavy hadrons, on the other hand, radial flow re- 
duces V2 at low Pt [10|- Assuming positive elliptic flow, 
v±{ip^O,TT) > v± (<p=f , the stronger transverse flow 
v± in the reaction plane pushes heavy particles to larger 
Pt more efficiently in the reaction plane than perpendic- 
ular to it. In extreme cases [lo| this can, for heavy parti- 
cles, even lead to a depletion of low-pT emission into the 
reaction plane when compared with out-of-plane emis- 
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FIG. 5: (Color online) Transverse momentum dependence 
of the elliptic flow parameter for pions and protons. Solid 
(dashed) lines are with (without) hadronic rescattering. 



sion, i.e. to a negative V2(j>t) at low pt (even though 
their pT-integrated total elliptic flow V2 is positive) . But 
even without going to extremes, this mechanism generi- 
cally reduces V2{pt) at low pT for heavy hadrons. So it 
is the generation of additional radial flow in the hadronic 
stage which is responsible for (most of) the mass-splitting 
of V2{pt) observed in the low pT region. 

This mechanism works even if the (extra) radial flow is 
not perfectly hydrodynamic, i.e. if (as is the case in the 
hadron cascade) the system does not remain fully ther- 
malized, with locally isotropic momentum distributions. 
Any type of anisotropic collective transverse motion will 
cause such a mass-splitting of 1)2 (pt) at low pT, as long 
as the hadron in question participates in the flow. It 
is worth mentioning that in hydrodynamic calculations 
about half of the final radial fiow in Au-|-Au collisions at 
RHIC is generated du ring the hadronic stage (see Fig. 7 
in [H^l and Fig. 5 in A similar increase in radial 

flow generated by the JAM cascade is documented in 
Fig. Ha). 

From these observations we conclude that the large 
magnitude of the integrated 1)2 and the strong mass or- 
dering of the differential V2{pt) observed at RHIC result 
from a subtle interplay between perfect fluid dynamics 
of the early QGP stage and dissipative dynamics of the 
late hadronic stage: The large magnitude of V2 is due to 
the large overall momentum anisotropy, generated pre- 
dominantly in the early QGP stage, whereas the strong 
mass-splitting between the slopes of V2{pt) at low pT 
reflects the redistribution of this momentum anisotropy 
among the different hadron species, driven by the contin- 
uing radial acceleration and cooling of the matter during 
the hadronic rescattering phase. 
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C. Spectra and elliptic flow for ^ mesons 



As noted in Sec. IIIDl </> mesons (consisting of strange 
quarks) have considerably smaller scattering cross sec- 
tions in JAM than non-strange hadrons (56j . They are 
therefore expected to show larger dissipative effects in 
our hybrid model and to not fully participate in the ad- 
ditional radial flow generated during the hadronic rescat- 
tering stage. In kinetic theory language, one expects that 
the (f> mesons decouple from rest of the system earlier 
than other, non-strange hadrons [131, thereby possibly 
opening a window to extract direct information on col- 
lective phenomena in the partonic stage from c/i-meson 
spectra fsoj . 

To study (j) mesons in our hybrid model we stabi- 
lize them by turning off their decay channels during the 
hadronic cascade. 
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FIG. 7: (Color online) Transverse momentum spectra for (j> 
mesons reconstructed from K^K~ decays in central (blue 
line), semi-central (red line) and peripheral (green line) 
Au+Au collisions, compared with PHENIX ^ and STAR 
[59I I data. Results from semi-central and peripheral collisions 
are divided by 10 and 100, respectively. Predictions from ideal 
hydrodynamics with Tth = 100 MeV are also shown as dashed 
lines. 



FIG. 6: (Color online) Normalized distribution of freezeout 
times for pions (dashed), protons (dotted), and (p mesons 
(solid) for |y| < 1 in Au+Au collisions at 6=2.0fm. 

Let us first check how early mesons decouple from 
the rest of the system. Figure |6] shows the normalized 
distribution of freeze-out times for pions, protons, and 
(j) mesons near midrapidity |y| < 1 in central collisions 
((fo)=2.0fm). Clearly, (j) mesons decouple earlier than 
pions and protons. The freeze-out time distribution for 
(j) mesons has a prominent peak at r=:8fm/c, roughly 
equal to the time of completion of QGP hadronization 
in hydrodynamic simulations. This indicates that only 
very few rescatterings happen for (j) mesons during the 
hadronic evolution. Similar results were obtained with 
the RQMD cascade in [s^l for ft baryons at SPS ener- 
gies and in [g^I for </> mesons and fl baryons at RHIC 
energies. The freeze-out time distributions for pions and 
protons are broadened by both clastic scatterings and res- 
onance decays. The long resonance decay tails of the dis- 
tributions arc important for interpreting the pion source 
function that was recently reconstructed by the PHENIX 



In Figure [71 pt spectra for (j) mesons from the hybrid 
model are compared with PHENIX [H and STAR [H 
data. Similar to the spectra for pions, kaons, and pro- 
tons in Fig.[3l we see good agreement with experiment at 
low Pt {pt < 1.5GeV/c). The discrepancy between our 
results and experiment at larger px may indicate the ap- 
pearance of a quark-antiquark recombination component 
in the intermediate pT region [63 . i63j . In the presence of 
such a component it is questionable to use the 0-meson 
spectra over the whole available pt region to extract the 
thermal freeze-out temperature and flow for (j) mesons 
[67t ; such a thermal model fit should be restricted to 
the region < 1-5 GeV/c even if data in that region are 
hard to obtain. 

In the hydrodynamic model simulations with Tth = 100 
McV, shown as dashed lines in Fig. [71 the 4> mesons pick 
up more additional radial flow during the hadronic stage, 
resulting in flatter p^-spectra than in the hybrid model 
and in the data in the low pt region. As we will show 
further below, better (j) data at lower pT and a simultane- 
ous analysis of the differential elliptic flow in this region 
should allow to further discriminate between different de- 
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scriptions of the hadronic rescattering stage. 

The effects of radial flow, and the difference in how ad- 
ditional radial flow generated during the hadronic rescat- 
tering stage is picked up by protons and mesons (which 
have rather similar masses), can be enhanced by study- 
ing the pt or transverse kinetic energy dependence of 
the (j)/p ratio. A thermalized medium without radial 
flow features mx-scaling, i.e. all mT^-spectra have iden- 
tical slopes, and for such a static flreball the (p/p ratio, 
when plotted as a function of transverse kinetic energy 
KEt = mT—mo, would be a constant horizontal line. For 
a thermalized expanding medium, mT-scaling is broken 
by radial flow (which couples differently to particles with 
different masses), resulting in a non-zero slope of the ra- 
tio (j)/p(KEiT). Perhaps somewhat counterintuitively, this 
slope of the (f>/p ratio does not grow monotonically with 
the radial flow v± but, after an initial rise, decreases 
again when the flow becomes so large that the hadron 
TOT-spectra become very flat; in the limit of "infinite 
flow" (i.e. 7j_ = l/y'l— oo) the hadron 7TiT-spectra, 
and thus their ratios, become again perfectly flat. 

In Figure [8] we show the (j>/p ratio, both as a function 
of transverse kinetic energy (right panel) and of pT (left 
panel) . It should be noted here that weak decay contribu- 
tion is not included in proton yields. In the latter case the 
connection to radial flow is less straightforward, since the 
kinematics of the transformation from niT to pt depends 
on mass and introduces additional growth with px for the 
ratio. In both representations one sees, however, by com- 
paring the curves for the hydro-|-cascade model without 
rescattering (corresponding to ideal hydrodynamics with 
Tth =169 MeV) and for the ideal hydrodynamic model 
with Tth = 100 MeV, that (i) the ratio increases with pr 
or KEt due to radial flow effects, and that (ii) the rate 
of increase drops when the freeze-out temperature Tth 



is decreased, due to build-up of additional radial flow. 
Surprisingly, the ratio increases even in pp collisions, but 
for entirely different reasons, unrelated to collective flow: 
The (j) spectrum from pp collisions shown in Fig. llOl below 
is considerably flatter than the proton spectrum, leading 
to the prominent rise of the (jy/p ratio with px. The most 
interesting feature of Fig.[8]is that the (p/p ratio from the 
hybrid model does not at all increase with px or KE^ (ex- 
cept at very low px < 500MeV/c). Instead, it decreases 
over almost the entire range of transverse kinetic energy 
shown in the figure. This decrease is due to the flat- 
tening of the proton spectrum by hadronically generated 
radial flow in which the weakly coupled (f> mesons do not 
participate. The comparison with pp collisions and hy- 
drodynamic model simulations in Fig. [5] shows that the 
observation of such a decreasing (p/p ratio would be an 
unambiguous signature for early decoupling of (p mesons 
from the hadronic rescattering dynamics. 

We now proceed to the discussion of dissipative effects 
during the hadronic rescattering stage on the differen- 
tial elliptic flow V2{px)- Figure [5] shows V2{px) from the 
hybrid model for tt, p, and p. We consider semi-central 
collisions (20-30% centrality), choosing impact parame- 
ter b = 7.2 fm. In the absence of hadronic rescattering 
we observe the hydrodynamically expected mass order- 
ing (px) > vl{px) > v^ipx) (Fig. IHa)), but just as in 
Fig. El (dashed lines) the mass splitting is small. Fig- 
ure [ni^b) shows the effects of hadronic rescattering: while 
the V2{px) curves for pions and protons separate as dis- 
cussed before (at low px the pion curve moves up while 
the proton curve moves down), V2{px) for the p meson re- 
mains almost unchanged [U . As a result of rescattering 
the proton elliptic flow ends up being smaller than that 
of the (p meson, (pr) < V2{px) for <pt < 1-2 GeV/c, 
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FIG. 9: (Color online) Transverse momentum dependence of the elliptic flow parameters for pions (dotted blue), protons 
(dashed green), and (j> mesons (solid red), for Au+Au collisions at &=7.2fm. (a) Before hadronic rescattering. (b) After 
hadronic rescattering. (c) Ideal hydrodynamics with Tth = 100 MeV. The results for pions and protons are the same as shown 
in Fig. [5] 



even though > nip. Hadronic dissipative effects are 
seen to be particle specific, depending on their scattering 
cross sections which couple them to the medium. The 
large cross section difference between the protons and (p 
mesons in the hadronic rescattering phase leads to a vi- 
olation of the hydrodynamic mass ordering at low pT in 
the final state. 



This is the most important new result of our work. 
Current experimental data (65l . [66| neither confirm nor 
contradict this predicted behavior, due to the difficulty 
of reconstructing low-pT mesons from their decay prod- 
ucts. If it turns out that high precision 0-meson V2 
data at low show violation of mass ordering, it will 
be evidence for strong momentum anisotropy having de- 
veloped already during the QGP stage, with the con- 
tribution carried by <f) mesons not being redistributed 
in PT by late hadronic rescattering. At intermediate 
Pt, recent data [H, [6^ confirm the prediction from the 
quark coalescence model [gl, [t^ that there the elliptic 
flow should scale with the number of constituent quarks: 
"^tiPT) ~ '^2'^ (pt) ~ |''^2(-Pt), in spite of the similar (p 
and p masses which are much larger than those of the pi- 
ons and K mesons. We hope that the present paper moti- 
vates an effort to extend these data to lower px in order to 
test our prediction here that, at low px, v^ipr) < vtiPT) 
in spite of > nip. While the former observation sug- 
gests that at intermediate pr {2GcV/c<pt <6GeV/c) 
quark coalescence during the quark-hadron phase tran- 
sition controls the finally observed elliptic flow of all 
hadrons, without measurable distortion by subsequent 
hadronic reinteractions, confirmation of our prediction 
would confirm the importance of hadronic rescattering on 
low-pT hadrons, with results that depend on the magni- 
tude of the scattering cross sections of the various hadron 
species. 



We close this section with a discussion of the implica- 
tions of our hybrid model results for the nuclear modifi- 



cation factor 



Raa{pi 



cLNaa 



coll prpcLpj-dy 



Taa 



(11) 



PTdpTdy 



The observed suppression of pion yields at intermediate 
to high Pt [i| provides evidence of jet quenching in rel- 
ativistic heavy-ion collisions. For baryons, this suppres- 
sion effect is counteracted in the intermediate pT region 
by collective flow effects which, at low pr, lead to a rise of 
the p/tt (or, more generally, heavy/light) ratio as a func- 
tion of Pt. Collective flow effects extend into the inter- 
mediate Pt region 2 GeV/c <pt < 6 GeV/c even though 
the hydrodynamic picture is known to gradually break 
down above px > 1.5 — 2.5 GeV/c [7l|. Quark coalescence 
is one of the key mechanisms by which \ow-pT collectiv- 
ity on the quark-gluon level is transferred to the hadron 
spectra at intermediate px during the hadronization pro- 
cess [62, 63, 69.], leading to (unsuppressed) values of Raa 
(or of i?cp , the ratio of yields per number of binary col- 
lisions in central and peripheral collisions) of order unity 
for baryons at pT ~ 2-3 GeV/c [l|, [6^. We will show 
that hadronic rescattering following QGP hadronization 
affects Raa at low vt instead. 

The PHENIX ^ and STAR [6^ Collaborations have 
recently measured i?cp for </) mesons. The PHENIX 
data show a suppression of 4> mesons by about a fac- 
tor 2 (with relatively large error bars) in the region 
1 GcV/c <pt < 3 GcV/c, consistent with that of pions, 
while protons and antiprotons are unsuppressed [5^. 
This seems to be in contradiction with collective flow 
arguments which predict i?^p > i?Qp since > nip, but 
consistent with the valence quark scaling predicted by 
the quark coalescence mode l |62| . [63l |69| . The more re- 
cent and precise STAR data [66[ , on the other hand, show 
an Rcp for (f) mesons that follows the one for pions and 
exceeds the one for protons for pT <l GeV/c, but then 
follows the rise of the proton Rcp above the pion one for 
PT>lGeV/c, lagging only slightly behind the protons 
and reaching a value halfway between pions and protons 
in the region pT ~ 2— 3 GeV/c where R^p peaks at a value 
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FIG. 10: (Color online) Invariant cross sections as a function 
of p t in non-singly diffractive pp collisions for pions, protons, 
[73 1 and (j> mesons. Dotted, dashed, and solid lines are results 
from PYTHIA for pions , protons, and (j) mesons, respectively. 



of - 1. 

Given this somewhat contradictory experimental situ- 
ation, we offer a prediction from our hybrid model (cau- 
tioning beforehand that this model does not include any 
quark-recombination contributions which are expected 
to become important above ^ 1-5— 2 GeV/c) in Fig- 
ure [TT] To construct this Figure, we first fitted the ex- 
perimentally measured pr-spectra for pions and protons 
[72| | as well as for (p mesons [73j in non-singly diffrac- 
tive (NSD) pp collisions {i.e. inelastic collisions exclud- 
ing single diffractive events). The fit, shown in Fig.lTOl is 
performed with the help of the event generator PYTHIA 
6.403 [zil which, once properly tuned, yields smooth ref- 
erence pT-spectra for pp collisions. PYTHIA is based on 
leading order perturbative QCD for semi-hard processes 
combined with a Lund string fragmentation scheme for 
soft particles. It works quite well f or p ions, protons and 
(j) mesons with default parameters [74| . except for a nec- 
essary readjustment of the K factor to K = 1.8. We note 
that here exceptionally this comparison includes all res- 
onance decays including weak ones since the STAR data 
show the inclusive spectra. We take the resulting spectra 
as our pp reference, after removing weak decay contribu- 
tions and multiplying them with the ratio crin/cNSD to 
correct for the NSD trigger. For the required cross sec- 
tions PYTHIA provides the estimates ctnsd = 32 mb and 
(Tin = 42 mb. 



With these reference spectra the nuclear modification 
factors Raa can now be calculated from the results 
shown in Figs. [H [3] and \7\ For pions, protons, and 
mesons they are shown as functions of pT in Fig. 1111 for 
Au-I-Au collisions at impact parameter 6 = 3.2fm {i.e. 0- 
10% centrality). Figure fTTKa) shows the predictions for 
the hybrid model. While for pions Raa{pt) is almost 
fiat, i?^^ — 0.15— 0.25, the Raa{pt) curves for protons 
and (f) mesons increase with pT as expected from ra- 
dial flow arguments (radial flow hardens the pt spec- 
tra for heavy particles). The rate of increase for the 
mesons is very similar to that for protons, culminating 
in a peak value of ~ 60% at pt ~ 1-2 — 1.4 GeV/c for 0's 
whereas the Raa for protons peaks at a value of — 60% 
near px^l-S GeV/c. FigurelTTTb) shows the correspond- 
ing curves for the ideal fluid dynamical simulation with 
Tth = 100 MeV. For pions and protons, the differences to 
the hybrid model are minor (at least in the pt range 
covered in the Figure) , reiterating the observation made 
in connection with Fig. [T] that the buildup of additional 
radial flow during the hadronic stage is similar in both 
models and viscous effects become clearly recognizable 
only at larger pT. For </> mesons one observes a much 
faster rise of RAAiPr) in the hydrodynamic approach, 
resulting in a larger peak value of ~ 105% at a larger 
Pt value (~ 1.7GeV/c) than for the hybrid model. The 
reason for these effects is obviously the larger amount of 
radial flow picked up during the hadronic stage in the 
hydrodynamic model and the resulting hardening of the 
(j) spectrum. The much weaker rise of R'^j^{pt) in the 
hybrid model can thus be traced directly to the lack of 
meson rescattering during the hadronic stage. 

We note that, even in the hydrodynamic model, the nu- 
clear modification factor Raa{pt) doesn't show a mono- 
tonic mass-ordering at low px. Naive expectations based 
on the mass-ordering of the spectral slopes (which re- 
flect radial flow effects) are invalidated by the fact that 
the (j) PT-spectra from pp collisions are flatter than the 
corresponding proton spectra. Since these spectra enter 
the denominator of Raa, they distort its pr dependence 
differently for protons and (f> mesons. 

We also comment that at ~ 2 GeV/c, the character- 
istics of the observed mass-scaling violation in Fig. [TlTa) 
are qualitatively similar to those expected (and observed) 
in the quark coalescence picture at intermediate pt 
(2 GeV/c <pT<6 GeV/c) . The differences are quan- 
titative: our prediction for Raa features neither a mono- 
tonic mass-ordering at low px nor the strict valence quark 
scaling predicted by the quark-coalescence picture at in- 
termediate Pt . 



IV. CONCLUSIONS 

We have studied effects of hadronic dissipation on the 
spectra, differential elliptic flow, and nuclear modifica- 
tion factor of pions, kaons, protons, and (j) mesons from 
Au-I-Au collisions at RHIC, using a hybrid model which 
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FIG. 11: (Color online) Nuclear modification factors Raa for pions (blue), mesons (red), and protons/antiprotons (green), 
for Au+Au collisions at & = 3.2 fm (corresponding to 0-10% centrality). Shown are predictions from (a) the hybrid model and 
(b) from ideal hydrodynamics with Tth = 100 MeV. 



treats the early QGP phase macroscopically as a per- 
fect fluid and the late hadronic phase microscopically 
with a hadronic cascade. For transverse momenta be- 
low 1.5GeV/c and not too peripheral collisions, the hy- 
brid model gives a reasonable description of the measured 
pion, kaon, proton and (f> meson pj^-spectra. In peripheral 
collisions (6 = 9 fm and larger) the model spectra tend 
to be somewhat steeper than measured. The centrality 
dependence of the differential elliptic flow V2{pt) of pi- 
ons, kaons and protons is better described by the hybrid 
model than in a purely hydrodynamic approach. 

For pions, kaons, and protons, which have relatively 
large scattering cross sections, hadronic rescattering is 
seen to generate additional collective transverse flow, but 
not so for the much more weakly interacting </> mesons. 
However, even for pions and protons the extra hadronic 
transverse flow effects are not "ideal" but exhibit ob- 
vious viscous features: Their p^-spectra are hardened 
while the growth of their elliptic flow V2{pt) with in- 
creasing pt is tempered by viscous corrections whose 
importance is in both cases observed to increase with 
transverse momentum. The well-known mass-splitting of 
the differential elliptic flow V2{pt) observed in hydrody- 
namic models is seen to be mostly generated during the 
hadronic rescattering phase and to be largely due to a 
redistribution of the momentum anisotropy built up dur- 
ing the QGP stage. This redistribution is caused by the 
mass-dependent flattening of the transverse momentum 
spectra by additional radial flow generated during the 
hadronic stage. The much more weakly interacting 
mesons do not participate in this additional radial flow 
and thus are not affected by this redistribution of momen- 
tum anisotropics: their differential elliptic flow remains 
almost unaffected by hadronic rescattering. The net re- 
sult of dissipative hadronic rescattering is therefore that 



the differential elliptic flow V2{pt) of protons drops below 
that of the (f> mesons, in violation of the hydrodynamic 
mass-ordering. A similar violation of the mass-ordering 
is seen in the nuclear modification factor Raa{pt) at 
Pt ^ 2 GeV/c where, after hadronic rescattering, the 
curve for (f) mesons ends up between those for pions and 
protons even though the </> is heavier than both of them. 
For the (j)/p ratio, the lack of interaction between the 
mesons and its accelerating hadronic environment should 
manifest itself in an unexpected but unambiguous de- 
crease with increasing transverse kinetic energy. 

The results presented here underscore the conclusion 
of Ref. [isl that hadronic dissipation may be very im- 
portant at RHIC and at lower beam energies and should 
be properly accounted for in attempts to quantitatively 
account for the experimental data collected from heavy- 
ion collisions. With V2{pt) and Raa(j>t) for \aw-pT 4> 
mesons and the dependence of the (ji/p ratio on px or 
transverse kinetic energy KE^, we have identified three 
additional critical observables which should be helpful in 
sorting out the interplay between hydrodynamic evolu- 
tion during the early QGP stage and dissipative hadronic 
expansion during the late stage of the hot and dense fire- 
balls created in these collisions. An accurate extraction 
of the value for the specific shear viscosity ry/s of the 
QGP created at RHIC requires a proper accounting for 
effects from late hadronic viscosity. Here, an attempt 
has been made to do this, by coupling the hydrodynamic 
model to a hadronic cascade. 
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